Fracturing highly disordered materials 
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We investigate the role of disorder on the fracturing process of heterogeneous materials by means 
of a two-dimensional fuse network model. Our results in the extreme disorder limit reveal that 
the backbone of the fracture at collapse, namely the subset of the largest fracture that effectively 
halts the global current, has a fractal dimension of 1.22 ± 0.01. This exponent value is compatible 
with the universality class of several other physical models, including optimal paths under strong 
disorder, disordered polymers, watersheds and optimal path cracks on uncorrelated substrates, hulls 
of explosive percolation clusters, and strands of invasion percolation fronts. Moreover, we find that 
the fractal dimension of the largest fracture under extreme disorder, df — 1.86 ±0.01, is outside the 
statistical error bar of standard percolation. This discrepancy is due to the appearance of trapped 
regions or cavities of all sizes that remain intact till the entire collapse of the fuse network, but are 
always accessible in the case of standard percolation. Finally, we quantify the role of disorder on the 
structure of the largest cluster, as well as on the backbone of the fracture, in terms of a distinctive 
transition from weak to strong disorder characterized by a new crossover exponent. 
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The brittle fracture of heterogeneous systems still rep- 
resents a major challenge from both scientific and techno- 
logical points of view. It has been the subject of intense 
scientific research in physics, material science, mechani- 
cal and metallurgical engineering The amount of 
small cracks and the shape of the largest fracture poten- 
tially determine whether or not the system still supports 
external loads and how catastrophic will be the rupture. 
Previous studies have shown that the degree of disorder 
in the material rules the transition between abrupt and 
gradual ruptures, i.e., the more heterogeneous is the ma- 
terial, more warnings one gets before the system collapses 
[B|. However, the nature of this transition, as well as the 
behavior of the fracturing system in the limit of strong 
disorder @, 0], remain as important open questions. Ev- 
idently, to model fracturing formation, it is fundamental 
to determine how stress is redistributed over the system 
while cracks appear, grow, and merge. Although concep- 
tually simple, the fuse network model is a good candidate 
for dealing with such problem, since it can clearly capture 
the essential features of the involved physical phenomena. 
In this model, resistors are used to mimic springs, in an 
approximate description for elasticity, where vector and 
tensor fields describing fracture mechanics are replaced 
by a scalar field representing the local strain [8 - 11 1 . 

The purpose of this Letter is to investigate the role 
of disorder on the scaling properties of the fuse network 
model at the critical collapse condition. We first study 
in detail the limiting case of extreme disorder. This 
is performed by measuring, through a purely geomet- 
rical technique, the fractal dimension of three different 
structures, namely the backbone of the fracture of broken 
bonds that halts current through the lattice, the largest 
fracture formed by the broken bonds (backbone plus dan- 



gling ends attached to it), and the total network of all 
broken bonds (largest cluster plus smaller clusters not 
attached to it). We then show how the self-similar be- 
havior of the resulting fracture topology crosses over from 
one regime to another as we move from weak to strong 
disorder. 

Let us start by describing the fuse network model 
@, EH- In the traditional version of this model, the elec- 
tric potential in a resistor network should provide a sim- 
plified description of the local strain in a fracturing sys- 
tem. A crack forms when the stress or, correspondingly, 
the electric current over a given resistor, surpasses a cer- 
tain threshold. Therefore, our system is a lattice where 
each bond is a resistor with a given conductance and fus- 
ing threshold value. For simplicity, we consider here the 
case in which the conductance is the same for all bonds, 
however, we expect similar results with a varying con- 
ductance We model the heterogeneity by assigning 
to each bond i a threshold given by tj = 10^^ , where Ri 
is a random number uniformly distributed in the interval 
— 1 < Ri < 1. Therefore, the distribution of thresholds is 
hyperbolic, P t {t) ~ t -1 , with upper and lower bounds 
given by 10^ and 10 _/3 , respectively. 

Generally speaking, we consider that a potential dif- 
ference is applied between two opposite sides of a re- 
sistor network, and solve Kirchhoff's law to determine 
the current passing through each bond 14|. Below the 



threshold value, each bond conducts according to Ohm's 
law, but once the current Ii at a given bond i reaches 
the threshold tj, the bond burns (or breaks, in the me- 
chanical terminology) and becomes an insulator. In this 
way, the largest current-threshold ratio, max(Ii/Ti), de- 
termines the next bond to be burnt. Here we assume that 
the potential starts from zero and raises slowly, allowing 
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FIG. 1: (Color online) The fuse network under extreme disor- 
der. In each step, a random bond belonging to the conducting 
part of the network is chosen to burn. Here the burnt bonds 
were removed and the red lines placed instead, but in the 
complementary lattice. The red lines represent cracks in the 
system and the green line close to the middle of the panel 
corresponds to a bond that will merge two large cracks, form- 
ing two cavities in the lattice (light blue and yellow dots). 
Since the current can not go through, none of the bonds in- 
side the cavities (light blue and yellow lines) can burn in the 
next steps. The dark blue and brown dots depict the limiting 
boundaries of each cavity. Note that the resulting fracture 
forms no loops. The purple dots in the middle top of the 
lattice correspond to sites connected to only one of the poles, 
therefore they do not participate in the transport, and the 
(purple) bonds connected to them do not burn either. 



only one bond to burn at each step. Pathological defects 
in the system are avoided by using a tilted square lat- 
tice, therefore in the first step all currents are the same, 
and the bond with the smaller threshold is the first to 
burn. In the following steps, inhomogencities are gradu- 
ally introduced in the system due to the burnt bonds, so 
that local threshold and current values should determine 
which one burns next. 

We initially focus on the case of a fuse network un- 
der extreme disorder, /3 — > oo, where thresholds are dis- 
tributed over many orders of magnitude. In this limit, for 
all practical purposes, one can assume that the smaller 
(larger than one) ratio between the thresholds of any two 
bonds in the lattice is larger than the largest ratio of the 
currents of any two bonds that constitute the conduct- 
ing part of the network. In this regime, any variability in 
the currents becomes irrelevant so that the next bond to 
burn should be the one with the smaller threshold among 
those that participate in transport. Since the thresholds 
are randomly distributed over the lattice, this is entirely 
equivalent to a process in which a random bond with fi- 
nite current, I > 0, is chosen to burn at each step. As 
the bonds burn, however, they may form a cavity that is 
connected to the rest of the lattice at a single node, as 
represented in Fig. [TJ Because all nodes inside a cavity 




FIG. 2: (Color online) Typical realizations of the fuse network 
at the point of disconnection (L = 32). In (a) f3 = 0.01 (weak 
disorder); in (b) ft = 1 (intermediate disorder); and in (c) 
j3 = 100 (strong disorder, for this lattice size). As represented 
in Fig. [1] the fractures occupy the' places of the burnt bonds, 
but in the complementary lattice. The black lines show the 
backbone with mass Mb, namely the number of burnt bonds 
forming the chain that effectively disconnects the bottom from 
the top of system. The set of all cracks connected to form the 
largest fracture, namely the backbone and its dangling ends 
(black and red lines, respectively) have mass Mf, while Mt 
denotes the total mass, i.e., the mass of the largest fracture 
plus all the other cracks (light blue lines) in the system. 



are equipotential, their connecting bonds are current-free 
and can not burn, regardless of their thresholds. 

In this way, the extreme disorder limit of the fuse net- 
work model can be viewed as a modified percolation prob- 
lem. In the standard percolation, bonds are randomly 
and sequentially occupied, while in the fuse network un- 
der extreme disorder, one only burns (occupies) bonds 
that belong to the conducting backbone of the lattice. 
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As shown in Fig. [TJ once a cavity is formed, the bonds 
inside remain unoccupied and, as a consequence, the clus- 
ters of occupied bonds in the complementary lattice do 
not form loops. Under this reasoning, the extreme dis- 
order limit of the fuse network model becomes a purely 
geometrical problem. Simulations performed with fuse 
networks at j3 = 10 8 confirm that these two models arc 
identical. Such a geometrical approach for the extreme 
disorder case greatly reduces the computational demand 
of the problem, therefore enabling us to simulate the frac- 
turing process for networks with linear size L going from 
16 to 2048, and using at least 1000 samples for the largest 
size. As shown in Fig. [21 we stop each realization when 
the two sides, where the potential difference is applied, 
become disconnected, i.e., no current can pass through 
the system. 

In Fig. [3] we show that the mass of the largest fracture, 
Mf, its backbone mass M b , as well as the set of all bro- 
ken bonds, M t , grow with the system size as power laws. 
First, the total mass of burnt bonds scales with the lin- 
ear size as, M t ~ L dt , with d t = 2.00 ± 0.01, suggesting 
that the bonds burn homogeneously through the lattice. 
The backbone grows as M b ~ L d » , with d b = 1.22 ± 0.01, 
which is statistically identical to the exponent obtained 
for optimal paths under strong disorder 15j, disordered 
polymers [1 61 ] . watersheds 11711 and optimal path cracks 
on uncorrelated substrates [l8| , the hulls of explosive per- 
colation clusters [19j . and the strands of invasion perco- 
lation fronts For the largest fracture, however, we 
obtain M f ~ L d f, with d f = 1.86 ± 0.01, which is dif- 
ferent from the fractal dimension of the largest cluster in 
2D percolation, d p = 1.8958 [H[22|. 
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In the limit of extreme disorder, the fracture backbone 
in the fuse model is identical to the one of loopless per- 
colation [23|. This can be seen by considering that the 
burning of fuses of a specific configuration of thresholds 
due to the extreme disorder just follows the sequence of 
their inverse rank, except if the fuse would close a loop or 
lie inside a nearly closed loop. In parallel, one can con- 
struct a configuration of loopless percolation by assigning 
an occupation rank to each site of the lattice identical to 
the ranking given by the thresholds (and obviously not 
occupying a site which would close a loop). The span- 
ning cluster of this percolation configuration then con- 
sists of the fracture of the fuse model and sites inside 
nearly closed loops that do not contribute to the back- 
bone. This shows that the bonds forming the backbone 
should be exactly the same in both models. It was previ- 
ously observed that the backbone of loopless percolation 
has the same fractal dimension as the optimum path in 
strong disorder [24[, therefore giving support to this ar- 
gument. 

The cavities can not change the fractal dimension of 
the backbone but may very well change the dimension of 
the largest fracture. The largest fracture comprises the 
backbone and the dangling ends attached to it. Once a 
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FIG. 3: (Color online) Size dependence results obtained for 
the fuse network model under extreme disorder. The total 
number of burnt bonds in the system scales as, Mt ~ L dt , 
with d t = 2.00 ± 0.01, indicating that the cracks appear ho- 
mogeneously in the 2D lattice. The largest fracture for this 
model grows as Mf ~ L df , with df = 1.86 ± 0.01, suggest- 
ing a different universality class for this structure than that 
of the spanning cluster in standard percolation [2lj]. Finally, 
the mass of the backbone of the largest fracture scales as 
M b ~ L db , with d b = 1.22 ± 0.01. The results shown in the 
inset confirm that the ratio between the masses of the largest 
clusters in the fuse and in the percolation models scales as, 



Mf/M p ~ L d f- d ", with d f - d p 



-0.03. 



cavity is formed in the largest fracture, it precludes the 
growing and attaching of other dangling ends inside it. 
As a consequence, the largest fracture in the extreme dis- 
order fuse model is certainly smaller than the largest clus- 
ter in standard percolation. Since we obtain, for the for- 
mer case, a fractal dimension smaller than the expected 
for standard percolation, we conclude that cavities form 
at every scale, and the difference between the two models 
becomes increasingly relevant as the system size grows. 

To further test our hypothesis, we performed simu- 
lations with a pinpoint algorithm that builds simulta- 
neously percolation clusters and fuse network fractures 
under extreme disorder. Precisely, at each step, we ran- 
domly choose a bond in the lattice and, if this bond is 
part of the conducting backbone, it is occupied in the 
percolation lattice as well as burnt in the fuse network. 
If this bond is inside one of the cavities, however, it is oc- 
cupied only in the percolation lattice and not in the fuse 
network. Once a cluster forms that crosses the system, 
the percolation condition is achieved in both models si- 
multaneously. If the fractal dimensions are in fact differ- 
ent, the ratio between the masses of the largest clusters, 
Mf/Mp, where M p is the mass of the spanning cluster in 
percolation, should vary with the system size as a power 
law, with an exponent given by df — d p w —0.03. The 
obtained results shown in the inset of Fig. [3] support our 
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FIG. 4: (Color online) The transition from weak to strong 
disorder. We show the mass of the largest fracture (a) and 
backbone (b) as function of the disorder parameter /3. In the 
week disorder regime, /3 < ft*, we have a linear growth for 
both backbone and largest fracture. As /3 increases the masses 
cross over for the strong disorder regime where a super-linear 
growth is observed. The exponents df and dt, controlling the 
scaling in the strong regime, are the same as in Fig $3$. The 
value fi* determining the onset of the crossover scales with 
the size of the system as f3* ~ L 9 . 



conjecture that the largest fracture of the fuse network 
under extreme disorder corresponds to a subset of the 
spanning cluster in percolation. Moreover, this new ob- 
ject is also self-similar, but with a slightly smaller fractal 
dimension. 

It is important to mention that some variants of the 
percolation model previously investigated also depart 
from the universality class of standard percolation. For 
example, in invasion percolation with trapping [251 ] , when 
a loop is formed in the growing cluster, no site or bond 
inside this loop can be occupied, constituting a trapped 
region. These trapped regions are somehow similar to our 
cavities, only that cavities are sections of the lattice out- 
side the conducting backbone of unoccupied bonds, while 
trapped regions are sections outside the infinite cluster of 
unoccupied bonds. In this case, statistically relevant de- 
viations have also be detected from the fractal dimension 
of the spanning cluster in standard percolation 26[ 27 1. 



Next we investigate how the behavior of the fuse net- 
work model crosses over from weak to strong disorder by 
gradually increasing the value of the parameter (3. Lo- 
cal currents are computed by applying Kirchhoff's law to 
each site of the network at each burning step, and solv- 
ing the resulting system of linear algebraic equations. In 
weak disorder, most of the burnt bonds belong to the 
backbone of the fracture that grows linearly with the 



system size, db 



1. As f3 increases, the fractal di- 



mensions obtained in the extreme disorder limit should 
be eventually recovered. Figure U shows how the masses 
Mf and Mb vary with the disorder parameter /3. For 
small values, (3 < 0.1, both masses are proportional to 
the system size, and depend weakly on ft. For intermedi- 
ate values, 0.1 < ft < /3*(L), the masses depend of (3 but 
still grow linearly with L, therefore indicating the persis- 
tence of the weak disorder regime. However, as one goes 
to larger values, (3 > f3*(L), the curves cross over to the 
strong disorder regime, where the masses show super- 
linear growth with system size, M t ~ L dt , Mf ~ L d f , 
and Mb ~ L db , but are again not dependent on (3. The 
value (3* marks the transition from weak to strong dis- 
order. One should expect that, above a characteristic 
length £, the system scales as in the weak disorder limit. 

Certainly, the length £ should depend on the strength of 

i 

the disorder, diverging in the extreme disorder as £ ~ (3 v . 
The onset of the strong disorder regime takes place when 
£ > L, resulting in f3* ~ L n . The insets of Figs. 0] show 
the data scaled by the fractal dimensions and f3* , with 
the controlling exponent 77 = 0.9 zb 0.1. The collapse of 
the curves in the transition region corroborates our anal- 
ysis and reveals the presence of a crossover between the 
two regimes. 

In conclusion, our results show that the threshold dis- 
order introduces a characteristic scale £ in the system. 
Below this scale, the fracture backbone displays a tor- 
tuous self-similar shape, with the same fractal dimen- 
sion of the optimum path under strong disorder 3, 13| 
and other previously investigated models [HI [l7l4l9l |. For 
scales larger than £, the fractures grow linearly with sys- 
tem size, consistent with the weak disorder regime. In 
the limit of extreme disorder, £ — > 00, the largest frac- 
ture has a fractal dimension of df = 1.86 ±0.01, close to, 
but different from the fractal dimension of percolation 
clusters (22)]. 
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